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Abstract 

The process 77 — > 4>J/ip is measured for <j)J/ip masses between threshold and 5 GeV/c 2 , us- 
ing a data sample of 825 fb _1 collected with the Belle detector. A narrow peak of 8.8I3' 2 , 
events, with a significance of 3.2 standard deviations including systematic uncertainty, is ob- 
served. The mass and natural width of the structure (named X(4350)) are measured to be 
(4350.6l£?(stat)±0.7(syst)) MeV/c 2 and (13^ 8 (stat) ± 4(syst)) MeV, respectively. The product 
of its two-photon decay width and branching fraction to <pJ/^p is (6.7jr 2 '4(stat) ± l.l(syst)) eV 
for J p = 0+, or (1.5l^(stat) ± 0.3(syst)) eV for J p = 2+. No signal for the Y(4140) -»■ (pJ/ip 
structure reported by the CDF Collaboration in B —> K + (f)J/if; decays is observed, and limits of 
r 77 (Y(414O))0(y(414O) -»■ <f)J/ip) < 41 eV for J p = 0+ or < 6.0 eV for J p = 2+ are determined 
at the 90% C.L. This disfavors the scenario in which the Y(4140) is a D* + D*~ molecule. 

PACS numbers: 14.40.Rt, 13.25.Gv, 13.66.Bc 
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In recent years, many new charmonia or charmonium-like states have been discovered. 
These states are not easily accommodated in the quark-model picture of hadrons In 
this Letter, we report the first investigation of the <pJ/ip system produced in the two-photon 
process 77 — > 4>J/ijj with the J/ip decaying into lepton pairs and <fi — > K + K~, to search for 
high mass states with J PC = ++ or 2 ++ , such as the tetraquark states and molecular states 
that are predicted by various models [2-4]. We find an unexpected structure in (j)J/ip mass 
near 4350 MeV/c 2 . 

In a related study of B + — > K + (pJ/ip decays, the CDF Collaboration reported evi- 
dence of a state called V(4140) with mass and width values of M = (4143.0 ± 2.9(stat) ± 
1.2(syst)) MeV/c 2 and T = (11.7l|jj(stat) ± 3.7(syst)) MeV Q. The Belle Collaboration 
searched for the V(4140) using the same B mode with a sample of 772 x 10 6 BB pairs (6j|. 
No significant signal was found although the upper limit on the production rate does not 
contradict the CDF measurement. 

There have been a number of different interpretations proposed for the V(4140), including 
a D* + D*~ molecule 0, [7-13], an exotic 1 _+ charmonium hybrid 0], a cess tetraquark 



state [3], or a natural consequence of the opening of the <f>J/i^ channel [14j. There are 
arguments against the interpretation of the K(4140) as a conventional charmonium state, 
such as the x'co or Xci [3; or a scalar D* + D*~ molecule since QCD sum rules [lfl, [l7[ 
predict masses that are inconsistent with the observed value. Assuming that the K(4140) 
is a D* + D*~ molecule with quantum numbers J PC = ++ or 2 ++ , Ref. predicts its 
two-photon width to be of order 1 keV, which can be tested experimentally at Belle. 

This analysis of 77 <fiJ/ip is based on a 825 fb _1 data sample collected with the Belle 
detector [l8| operating at the KEKB asymmetric-energy e + e~ collider (l9l ] . About 90% of 
the data were collected at the T(nS) (n = 1,3,4,5) resonances, and about 10% were taken 
at a center-of-mass (CM.) energy that is 60 MeV below the T(45) peak. 

The detector is described in detail elsewhere [llj]. It is a large-solid-angle magnetic 
spectrometer that consists of a silicon vertex detector (SVD), a 50-layer central drift chamber 
(CDC), an array of aerogel threshold Cherenkov counters (ACC), a barrel- like arrangement 
of time-of-flight scintillation counters (TOF), and an electromagnetic calorimeter comprised 
of CsI(Tl) crystals (ECL) located inside a superconducting solenoid coil that provides a 1.5 T 
magnetic field. An iron flux-return located outside of the coil is instrumented to detect K\ 
mesons and to identify muons (KLM). 

We use the program treps [20] to generate signal Monte Carlo (MC) events. In this gen- 
erator, the two-photon luminosity function is calculated and simulated events are generated 
at a specified fixed 77 CM. energy (VK 77 ) using the equivalent photon approximation (2lj . 
The efficiency for detecting 77 — » X — > <pJ/ip — > K + K~£ + £~ (£ = e, /z) is determined by 
assuming J p = + or 2 + and a zero intrinsic width for the X. 

We require four reconstructed charged tracks with zero net charge. For these tracks, 
the impact parameters perpendicular to and along the beam direction with respect to the 
interaction point are required to be less than 0.5 and 4 cm, respectively, and the transverse 
momentum in the laboratory frame is restricted to be higher than 0.1 GeV/c. For each 
charged track, information from different detector subsystems is combined to form a like- 



lihood Ci for each particle species 22] . Tracks with TZk = r fr > 0-6, are identified as 
kaons with an efficiency of about 97% for the tracks of interest; about 0.4% are misidentified 
tracks 22|. For electron identification, the likelihood ratio is defined as TZ e 



7T 



where C e and C x are the likelihoods for electron and non-electron, respectively, determined 
using the ratio of the energy deposit in the ECL to the momentum measured in the SVD 
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and CDC, the shower shape in the ECL, the matching between the position of charged 
track trajectory and the cluster position jn the ECL, the hit information from the ACC 

For muon identification, the likelihood ratio 
and Ck are the likelihoods for muon, pion 



and the dE/dx information in the CDC [23 
is defined as 1Z, 



where C, 



and kaon hypotheses, respectively, based on the matching quality and penetration depth of 
associated hits in the KLM (HI]. For electrons (muons) from J/ip decay, both of the tracks 
should have TZeilZ^) > 0.1. The lepton ID efficiency is about 99% for J/ip — > e + e~ and 90% 
for J /if; — > There are a few background events due to photon conversions with the 

conversion leptons misidentified as kaon candidates in the e + e~ mode; these are removed by 
requiring lZ e < 0.75 for the kaon candidates. 

The magnitude of the vector sum of the four tracks' transverse momenta in the CM. 
frame, which approximates the transverse momentum of the two-photon-collision 

system, is required to be less than 0.2 GeV/c in order to reduce backgrounds from non-two- 
photon processes and two-photon-processes with extra particles other than <ft and J/ip in 
the final states. Figure [1] shows the | 52 P% | distributions from the data and a MC signal 
simulation before the | X)P t *| requirement. Here the K + K~ and £ + £~ invariant masses are 
required to be within the </> and J / ip signal regions, respectively. 
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FIG. 1: The magnitude of the vector sum of (fiJ/ijj transverse momenta with respect to the beam 
direction in the e + e~ CM. frame for the selected (fiJ/ijj events. Points with error bars are data. 
The dot-dashed, solid, and dotted histograms are MC simulations for 77 — > (pj ' /ip with the (pJ/ijj 
mass fixed at 4.20, 4.35, and 4.50 GeV/c 2 , respectively (normalized to the number of events with 
I P* I < 0.2 GeV/c). The arrow shows the position of the | Y^P*\ requirement. 



A scatter plot of M(£ + £~) versus M(K + K~) for the selected K + K~£ + £~ events is shown 
in Fig. [21 where we can see clear J/ip and <fi signals. A partial correction for final state 
radiation and bremsstrahlung energy loss is performed by including the four-momentum of 
every photon detected within a 50 mrad cone around the electron and positron direction 
in the e + e~ invariant mass calculation. We define a J/ip signal region as 3.077 GeV/c 2 < 
m e + e - < 3.117 GeV/c 2 (the mass resolution is about 10 MeV/c 2 ), and J/ip mass sidebands 
as 3.0 GeV/c 2 < m £+£ - < 3.06 GeV/c 2 or 3.14 GeV/c 2 < m e+e - < 3.20 GeV/c 2 . We also 
define a <fi signal region as 1.01 GeV/c 2 < mx+K- < 1-03 GeV/c 2 (the full width at half 
maximum (FWHM) of the <fi signal is 5.9 MeV/c 2 ), and 4> mass sidebands as 1.00 GeV/c 2 < 
m K+K - < 1.01 GeV/c 2 or 1.03 GeV/c 2 < m K+K - < 1.08 GeV/c 2 . 
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FIG. 2: A scatter plot of M{l + i~) versus M(K + K~) for the selected K + K~£^ 
size of the boxes is proportional to the number of events. 



events. The 



Figure |3] shows the (pJ/ip invariant mass distribution [25|], together with the background 
estimated from the normalized J/ip and <p mass sidebands. No K(4140) signal is evident. 
Assuming that there is no background within the K(4140) mass region and the number of 
signal events follows a Poisson distribution with a uniform prior probability density function, 
a Bayesian upper limit on the number of the K(4140) signal events is estimated to be 2.3 



at the 90% CL. |26j. However, there is a clear enhancement at 4.35 GeV/c 2 , where the 
background level estimated from the normalized J/ip and <p mass sidebands is very low. 
Other possible backgrounds that are not included in the sidebands, such as 77 — > <pJ/ip + X 
and e + e~ — > <pJ/ip + X where X may indicate one or more particles, and 77 — > <fiJ /ip with 
the J/ip and <p decaying into final states other than lepton pairs and K + K~, are found to 
be very small after applying all of the event selection criteria. 
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FIG. 3: The (ft J /if) invariant mass distribution of the final candidate events. The open histogram 
shows the experimental data. The fit to the 4>J/ip invariant mass distribution from 4.2 to 5.0 
GeV/c 2 is described in the text. The solid curve is the best fit, the dashed curve is the background, 
and the shaded histogram is from normalized (j) and J/ip mass sidebands. The arrow shows the 
expected position of the y(4140). 



In order to obtain resonance parameters for the structure at 4.35 GeV/c 2 , an unbinned 
extended maximum likelihood method is applied to the <pJ/ip mass spectrum in Fig. [31 The 
distribution is fitted in the range 4.2 to 5.0 GeV/c 2 with an acceptance-corrected Breit- 
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Wigner (BW) function convoluted with a double Gaussian resolution function as the signal 
shape and a constant term as the background shape. The shape of the double Gaussian 
resolution function is obtained from MC simulation. 

From the fit, we obtain 8.8^2 signal events, with a mass m = (4350. 6i 5 j) MeV/ c 2 and a 
width T = (13lg 8 ) MeV. The statistical significance of this structure is estimated to be 3.9a, 
from the difference of the logarithmic likelihoods, —2 ln(L /L max ) = 21, taking the difference 
in the number of degrees of freedom (Andf = 3) in the fits into account, where L and L max 
are the likelihoods of the fits with and without a resonance component, respectively. In the 
following we refer to it as the X(4350). The significance of the signal decreases to 3.2cr if a 
linear function is used to model the background shape in the fit. 

We use an ensemble of simulated experiments to estimate the probability that background 
fluctuations alone would produce signals as significant as that seen in the data. We generate 
4>J/ijj mass spectra based on a uniform distribution alone with 24 events, the same as 
observed in data, and search for the most significant fluctuation in each spectrum in the 
mass range from 4.2 to 5.0 GeV/c 2 , with widths in a range between 3 MeV (half of the 
resolution) and 130 MeV (ten times the observed width). From these spectra we obtain 
the distribution for — 2 ln(L /I/ ma x) in pure background samples, and compare it with the 
signal in the data. We performed a total of 0.5 million simulations and found 65 trials 
with a — 2 ln(Lo/I/ max ) value greater than or equal to the value obtained in the data. The 
resulting p value is 1.3 x 10 -4 , corresponding to a significance of 3.8a. Generating events in 
a wider 4>J/ip mass range, or fitting with different width range would change the resulting 
significance, but the dependence is weak for a signal as narrow as the X(4350). 

The product of the two-photon decay width and branching fraction is obtained using 
the formula: T 77 (i?)i3(-R — > final state) = N/[(2J + l)e/C£ int ], where N is the number of 
observed events, e is the efficiency, J is the spin of the resonance, and £j nt is the integrated 
luminosity. K, is a factor that is calculated from the two-photon luminosity function £ 77 (M^) 
for a resonance with mass Mr using the relation: K, = 47t 2 £ 77 (Mr)/M^, which is valid when 
the resonance width is small compared to its mass (widths are smaller than 1% of the masses 
in the V(4140) and V(4350) cases). The /C parameter is calculated to be 0.46 fb/eV and 
0.36 fb/eV for the V(4140) and the V(4350), respectively, using treps Q. The efficiencies 
are 0.30% and 0.41% for J p = 0+ and 2+ respectively, at 4.143 GeV/c 2 , and 7.90% and 
6.98% for J p = + and 2 + respectively, at 4.35 GeV/c 2 . From the above values, we obtain 
r 77 (y(4140))£(V(4140) (f)J/ip) < 36 eV for J p = 0+ or < 5.3 eV for J p = 2+ at the 
90% C.L., andr 77 (X(4350))23(X(4350) -> <j>J/if>) = (6.7±|;1) eVfor J p = 0+ or (1.5i^) eV 
for J p = 2 + , where the errors are statistical only. 

There are several sources of systematic errors for the measurements of the products of the 
two-photon decay width and branching fractions. The particle identification uncertainties 
are 1.2%/kaon and 0.8%/lepton. The uncertainty in the tracking efficiency for tracks from 
J/i/j decays is 1% per track, while that for kaon tracks from decays range from 2.4% to 
1% per track as the average transverse momentum increases from 0.15 to 0.3 GeV/c. The 
efficiency uncertainties associated with the J/ip and <fi mass requirements are determined 
from the studies of the very pure e + e~ —tif)'—t 7i + 7i^J/ip [27| and e + e~ — > 07r + 7r~ [28| event 
samples. The detection efficiencies for J ftp and <fi mesons are lower than those inferred from 
the MC simulations by (2.5 ± 0.4)% and (2.0 ± 0.5)% relatively, respectively. We take 
0.96 as the efficiency correction factor, and 0.7% is included in the systematic error due 
to the J/tp and <fi mass requirements. The statistical errors in the MC samples are 2.3% 
and 0.9% for V(4140) — > <pJ/ip and V(4350) — > <f>J/ip, respectively. The accuracy of the 
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two-photon luminosity function calculated by the treps generator is estimated to be about 
5% including the error from neglecting radiative corrections (2%), the uncertainty from the 
form factor effect (2%), and the uncertainty in the total integrated luminosity (1.4%) [20]. 
The trigger efficiency for four charged track events is rather high because of the redundancy 
of the Belle first level multi-track trigger. According to the MC simulation, the trigger and 
preselection efficiency for the final state has little dependence on the 4>J/ip invariant mass, 
with an uncertainty that is smaller than 5%. From Ref. [26f, the uncertainty in the world 
average values for B((f> -> K+K~) is 1.2% and that for B(J/if) £ + £~) = B(J/ip e+e~) + 
B(J/tp is 1% where we have added the errors of the e + e~ and modes linearly. 

The uncertainty in the yield of X(4350) signal events due to the <f)J/ip mass spectrum fit is 
estimated to be 15% by varying: the order of the background polynomial (14%), resonance 
parameterization (0.5%) and the (frJ/ip mass resolution (0.9%). Assuming that all the sources 
are independent and adding all uncertainties in quadrature, we obtain the total systematic 
errors on r 77 (V(4140))£(V(4140) <j>J/ip) and r 77 (X(4350))fi(X(4350) -»■ (f>J/ip) to be 
12% and 17%, respectively. 

For the systematic errors in the X(4350) mass and width, the uncertainties in the mass 
resolution (0.1 MeV/c 2 and 0.9 MeV), the parameterization of the resonance (0 MeV/c 2 
and 0.3 MeV) and the background shape (0.7 MeV/c 2 and 3.9 MeV) are considered. As- 
suming that all the sources are independent and adding them in quadrature, we obtain the 
total systematic errors on the X(4350) mass and width to be 0.7 MeV/c 2 and 4.1 MeV, 
respectively. 

In summary, we report results of the first search for V(4140) — > 4>J/ip in the two-photon 
process 77 — > <f)J/ip. No V(4140) signal is observed, and upper limits on the product of the 
two-photon decay width and branching fraction of Y(4140) — > <f>J/ip are established to be 
r 77 (V(4140))£(V(4140) -> (f)J/ip) < 41 eV for J p = 0+ or < 6.0 eV for J p = 2+ at the 
90% C.L. In the determination of the r 77 (V(4140))i3(V(4140) <pJ/tp) upper limits, the 
efficiencies have been lowered by a factor of 1 — u sys to obtain a conservative estimate, where 
o sys is the total relative systematic error. The upper limit on r 77 (V(4140))i3(V(4140) — > 
(ftJ/i/j) from this experiment is lower than the prediction of (176lg3 7 ) e V for J PC = ++ , 
(1891^) eV for J PC = 2++ (calculated by us using the values in Ref. [2| and total width of 
the V(4140) from CDF Q). This disfavors the scenario in which the V(4140) is a D* + D* 
molecule with J PC = ++ or 2 ++ . 

We find evidence for an unexpected new narrow structure at 4.35 GeV/c 2 in the 
4>J/ip mass spectrum with a significance of 3.2 standard deviations including system- 
atic uncertainty. If this structure is interpreted as a resonance, its mass and width are 
(4350.6iti(stat) ± 0.7(syst)) MeV/c 2 and (13i^ 8 (stat) ± 4(syst)) MeV, respectively. The 
product of its two-photon decay width and branching fraction to <pJ/ip is measured to 
be r 77 (X(4350))£(X(4350) ->• <j>J/if>) = (6.7l|l(stat) ± l.l(syst)) eV for J p = 0+ or 
(1.5i[J;g(stat) ±0.3(syst)) eV for J p = 2 + . We note that the mass of this structure is consis- 
tent with the predicted values of a cess tetraquark state with J PC = 2 ++ in Ref. jl] and a 
D* + D*q molecular state in Ref. Q. In a recent paper 29 1, the possibility that the A (4350) 
could be an excited P-wave charmonium state, x'c2> was a l so discussed. 

We thank the KEKB group for excellent operation of the accelerator, the KEK cryogen- 
ics group for efficient solenoid operations, and the KEK computer group and the Nil for 
valuable computing and SINET3 network support. We acknowledge support from MEXT, 
JSPS and Nagoya's TLPRC (Japan); ARC and DIISR (Australia); NSFC (China); DST (In- 
dia); MEST, KOSEF, KRF (Korea); MNiSW (Poland); MES and RFAAE (Russia); ARRS 



(Slovenia); SNSF (Switzerland); NSC and MOE (Taiwan); and DOE (USA). 



[1] See for example: E. S. Swanson, Phys. Rept. 429, 243 (2006); E. Klempt and A. Zaitsev, 
Phys. Rept. 454, 1 (2007); S. Godfrey and S. L. Olsen, Ann. Rev. Nucl. Part. Sci. 58, 51 
(2008). 

[2] T. Branz, T. Gutsche, and V. E. Lyubovitskij , Phys. Rev. D 80, 054019 (2009). 

[3] Fl. Stancu, larXiv:0906.2485l 

[4] J. R. Zhang and M. Q. Huang. larXiv:0905.4672l 

[5] T. Aaltonen et al (CDF Collaboration), Phys. Rev. Lett. 102, 242002 (2009). 

[6] C. Z. Yuan (BES and Belle Collaborations). larXiv:0910.31381 

[7] X. Liu and S. L. Zhu, Phys. Rev. D 80, 017502 (2009). 

[8] G. J. Ding, Eur. Phys. J. C 64, 297 (2009). 

[9] N. Mahajan, Phys. Lett. B 679, 228 (2009). 
[10] X. Liu and H. W. Ke, Phys. Rev. D 80, 034009 (2009). 
[11] J. R. Zhang and M. Q. Huang. [arXiv:0905.4178[ 

[12] R. M. Albuquerque, M. E. Bracco and M. Nielsen, Phys. Lett. B 678, 186 (2009). 

[13] R. Molina and E. Oset, larXiv: 0907.30431 

[14] E. van Beveren and G. Rudd. larXiv: 0906 .22781 

[15] X. Liu, Phys. Lett. B 680, 137 (2009). 

[16] Z. G. Wang, Eur. Phys. J. C 63, 115 (2009). 

[17] Z. G. Wang, Z. C. Liu and X. H. Zhang, Eur. Phys. J. C 64, 373 (2009). 
[18] A. Abashian et al. (Belle Collaboration), Nucl. Instr. and Methods Phys. Res. Sect. A 479, 
117 (2002). 

[19] S. Kurokawa and E. Kikutani, Nucl. Instr. and Methods Phys. Res. Sect. A 499, 1 (2003), 

and other papers included in this volume. 
[20] S. Uehara, KEK Report 96-11 (1996). 
[21] C. Berger and W. Wagner, Phys. Rept. 146, 1 (1987). 
[22] E. Nakano, Nucl. Instr. and Meth. Phys. Res. Sect. A 494, 402 (2002). 
[23] K. Hanagaki et al, Nucl. Instr. Meth. A 485, 490 (2002). 
[24] A. Abashian et al, Nucl. Instr. Meth. A 491, 69 (2002). 

[25] In this Letter, m K + K -£+£- — mg+£- — m K + K - + mj/^ + is used instead of the invariant 
mass of the four final-state particles to improve the mass resolution. Here rrij^ and are 
the nominal masses of J/ip and <f> [2a ]. 

[26] C. Amsler et al. (Particle Data Group), Phys. Lett. B 667, 1 (2008) [Page 328, Table 32.3]. 

[27] C. Z. Yuan et al. (Belle Collaboration), Phys. Rev. Lett. 99, 182004 (2007). 

[28] C. P. Shen et al. (Belle Collaboration), Phys. Rev. D 80, 031101(R) (2009). 

[29] X. Liu, Z. G. Luo and Z. F. Sun, larXiv:0911.3694l 



9 



